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Abstract 

Two possible interpretations of FRW cosmologies (perfect fluid or 
dissipative fluid) are considered as consecutive phases of the system. 
Necessary conditions are found, for the transition from perfect fluid to 
dissipative regime to occur, bringing out the conspicuous role played 
by a particular state of the system (the "critical point"). 

1 Introduction 

It is already a well stablished fact that a variety of line elements may sat- 
isfy the Einstein equations for different (physically meaningful) stress-energy 
tensors (see [|1| and references therein). 

Particularly interesting is the situation with Friedmann- Robertson- Walker 
(FRW) models, which, as it has been shown and references therein), 

do not necessarily represent perfect fluid solutions, but can also be exact 
solutions for viscous and heat conducting fluids, with or without an electro- 
magnetic field. 

Thus, for example, it has been shown 0, that the zero- curvature FRW 
metric 

ds 2 = -dt 2 + R 2 (t) (dr 2 + r 2 d6 2 + r 2 sin 2 Odtf) (1) 
satisfies Einstein equations, not only for perfect fluid matter distributions 

T,u, = (p + p) V„ + pg^, (2) 

but also, for the stress-energy tensor of a magnetohydrodynamic viscous fluid 
with heat conducion, viz. 

Tp, = E tlv + (p + p- i®)u^u u + (p - £6)5^ - 2r]a^ + q^Uv + q v u^ (3) 

where E^, p, p, £, rj, O, and g M denote respectively, the electromagnetic 
stress-energy tensor, the energy density, the pressure, the bulk and shear 
viscosity coefficients, the expansion, the shear tensor and the heat flux vector. 

A fundamental difference between both interpretations of the models re- 
sides in the fact that observers moving along (perfect fluid case) are comov- 
ing with the fluid, whereas those moving along w M (dissipative case) detect 
a radial velocity of fluid particles (tilted models) Q. In other words, in the 
dissipative models the source is necessarily non-comoving or tilting relative 
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to the geometrically "preferred" Ricci eigenvector [31. According to Coley 
and Tupper this tilting velocity might be related to the observed motion 
of our galaxy relative to the microwave background radiation (CMB), and 
seems to be consistent with measurements of the dipole anisotropy of the 
CMB temperature distribution ||. 

A detailed analysis of many of those models shows that their observational 
predictions are in good agreement with available data ||. 

In this work, both interpretations are not considered as alternative, mutu- 
ally exclusive possibilities, but rather as consecutive phases in the evolution 
of the model. More specifically, we shall consider the case when, initially, the 
system evolves as perfect fluid without dissipation, and then at some instant 
starts to dissipate, tending eventually again to a perfect fluid as t — > oo. The 
question we want to answer here is: what could be the necessary condition 
for that transition to occur?. 

For simplicity we shall consider neither electromagnetic field {E^ u = 0) 
nor bulk viscosity (£ = 0). Furthermore, the viscosity term will be modeled 
through an anisotropic fluid and only transport equation for the heat flux 
will be needed. The reason of this simplification becomes clear later. Also, it 
should be kept in mind that the approach used in this work only provides 
a "snapshot" of the system at the moment it enters the dissipative regime, 
thereby giving only a hint about its tendency, and not a complete description 
of its evolution. All these simplifications represent the price to pay in order 
to deal with mathematically tractable equations. 

Two different initial situations will be considered: 

1. The system initially evolves strictly as a perfect fluid. By this we 
mean that there is not dissipation whatsoever and preferred observers 
(associated with the four velocity of the fluid) are (strictly) comoving 
with the fluid. 

2. The system initially evolves very close to the perfect fluid regime, but 
there is still an small (in a sense to be defined below) amount of dissi- 
pation, because of which, preferred observers are not strictly comoving 
with the fluid, but detect an small radial motion of the fluid, so that 
quadratic and higher powers, as well as time derivatives, of radial ve- 
locity are neglected during this "quasi-perfect fluid" regime. 

Then it is assumed that at some instant (say t = t), our system starts to 
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deviate from either of the initial regimes and enters into dissipative regime. 
In the first case the system leaves the perfect fluid condition by starting to 
dissipate. In the second case the system abandons the quasi-perfect fluid 
regime by increasing the rate of dissipation. 

We shall find conditions for transitions described above to occur. With 
this purpose we shall evaluate field equations and transport equations im- 
mediately after t, where immediately means on a time scale of the order 
of relaxation time. We shall see that the so called "critical point" plays a 
fundamental role in the occurrence of the transition to the dissipative regime. 

The critical point corresponds to the situation when (in relativistic 

units) 



where k, T and r denote heat conduction coefficient, temperature and relax- 
ation time respectively. 

It has been shown || that under that condition, the effective inertial 
mass of any fluid element, just after the system departs from equilibrium, 
vanishes. Also, it has been shown that as a self-gravitating system approaches 
the critical point, its active gravitational mass behaves in such a way that 
enhances the instability of the system ||. 

It will be shown below that transition from perfect fluid to dissipative 
regime is allowed only if the system is at the critical point. However, transi- 
tion from the quasi-perfect fluid phase to a dissipative regime is, in principle, 
always possible, and closer is the system to the critical point, faster is the 
transition. 

2 Field and Transport equations 
2.1 Field equations 

First of all, note that for any locally-Minkowskian observer, comoving with 
the fluid, the four velocity takes the obvious form 




(5) 
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then, performing a Lorentz boost to the frame with respect to which a fluid 
element has radial velocity u>, we find (in the coordinate system of ([!])) 

u»=(^J==,^EL,0,o) (6) 

\VT^ J 

which is the same expression (2.1) of 0, with their a's and /3's defined in 
terms of u by 



Next, the heat flux vector g M , satisfying q^u^ = 0, is given by 

g M = gf^=,^£=, 0,0^ (9) 



which is the same expression (2.3) of 0, with Q 2 = 57^- 

Then the energy-momentun tensor (^j of the source may be writen as 

T^ u = (p + P±)u IM u u + P±g fMU + (P r - Pl)s m s„ + QfjUu + (10) 

where P r and Pj_ denote the two principal (anisotropic) stresses, linked to p 
and r]a^ u by relations specified below, and the vector s M satisfies conditions 

s% = l ; s% = ; T^ l s v = P r (11) 

Then, Einstein equations read (in relativistic units and omitting the 8tt 
factor) 

- w - = T -^lp + P T u 2 -2Qu] (12) 
= (p + P r )u - Q(l + uj 2 ) (15) 
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where dots denote derivative with respect to t. 

Comparing (Q) and ([H]) the following relations follow 

P r =p-^a lx {l-u 2 ) (16) 
^=P+§^ii(l-^ 2 ) (17) 

2 U) 2 v _U 2 CTn 2U 2 (T 2 2 UJ 

a °° - 3 JT^) " " ~r 2 R 2 (l-u 2 ) ~ ~R ao1 (18) 

where X is given by 

X " (1-^)3/2 ( 19 ) 

and prime denotes derivative with respect to r. 
Observe that from (|i"3|) and (|i4|) , one obtains 

(p + P ± V 2 + (P r - P ± ) -2Quo = (20) 

and using (p~5f ) and (p0| ) it follows 

(Pr " Pi 



+ g(i-w 2 ) = o (2i) 



Thus, if P r = Pj_ (which implies by virtue of (|!6|)-(|l8l), rja^ = 0), then 
to — 1, which is obviously inadmissible. Therefore, pure heat conduction 
(plus a fluid) is not consistent with FRW metric. 

After a simple algebra, field equations (|T2]) — dl5|) may be written as 

P=^f + (Pr-P±) (22) 
(2R R 2 \ 

P± = ~ IT + 155 ( 23 ) 



P r = -(f + §)+g. (24) 
g (1-^=2^(1-1) (25) 
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It is worth noticing that (p3"|)-(f2"5|) imply that dissipation takes place only if 
u ^ and P r ^ P±. 

Later it will be convenient to write R(t) in the form 

R(t) = t n (26) 

where in general n = n(t). Thus, n = 2/3 reproduces the Einstein-de Sitter 
model, whereas 

n = - — (27) 
hit v ' 

reproduces the exponential inflationary model with R ~ e Hot , p + p = and 

P = Pr = P±- 

Using (|26|), field equations (|22|)-([25|) become 

lnt 3n 2 

p = 3 (hint) 2 + Qnh— + — + (P r - P±) (28) 

, , , 2 hit 3n 2 , 4h 2n , . 

P ± = -3 (nlnt) - Qnh - - 2nlnt + — (29) 

- y J t t 2 t t 2 y 1 

P r = Quo-3 (hint) 2 - Qnh— - ^- - 2nlnt - — + ^ (30) 
Q(l-uj 2 )=uj(^-^-2hlnt ] j (31) 

2.2 Transport equations 

In order to avoid the important drawbacks of Eckart's theory, we shall use 
transport equations derived from the Israel- Stewart approach 0- 11|. How- 



ever, since these equations are going to be evaluated immediately after the 
system leaves the perfect fluid (or the quasi-perfect fluid) regime, then the 
truncated version of the theory leads to the same result as the full theory, 
and therefore the former will be used here. For the same reason, the shear 
viscosity contributions result in terms which are negligible in the approxima- 
tion used here (see below). Therefore we shall deal only with the transport 



equation for the heat flux, which reads |TT 



<f + Th^u<qp. a = — k (h^Tp + Ta M ) (32) 
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where as usual, denotes the projector onto the three space orthogonal to 
w M , and a M denotes the four acceleration. 

In terms of n, the only non vanishing independent component of (^) is 
the r-component, which reads 



Q 



«{ 



+ T 



t n {l-uj 2 ) 1/2 
ujf 



+ r 



+ 



Q 



Q'uj 



t n (l-LU 2 ) 

V 



t 2n (1 - UJ 2 



t n (l-u 2 ) t 2n (l-u 2 ) 



UJ (jJOJ 

-J 2 f ~ " ' ^ 



t 2n n _ ^2} 



+ 



UJ 



t n a - uo 2 



n In t + — 



} (33) 



where prime denotes derivative with respect to r (the t-component of fl32|) is 
just fl33|) multiplied by to). 



3 Entering into dissipative regime 

As was already mentioned, an essential feature of the dissipative regime, 
as implied by ( P5| ) or ([31]), is the non comoving nature of observers of the 
congruence m m (w / 0). This means that u may be used as control variable 
to assess how far (or close) is the system from the perfect fluid regime. 

So, let us now assume that for t < t, our system is in either the perfect 
fluid regime or in the quasi-perfect fluid regime. Both of which are charac- 
terized as follows: 

1. Perfect fluid regime 

Q = lu = (34) 
Preferred observers are strictly comoving with the fluid. 

2. Quasi-perfect fluid regime 

oo^ 0(e) (35) 
Q w u w to 2 w Quj w g 2 w (36) 

Preferred observers are not strictly comoving with the fluid, but (|35|) 
and ( |36"D hold, i.e. within this regime the system always remains "close" 
to the perfect fluid phase. In both cases, of course, r derivatives of Q 
and uj, are of the same order of magnitude as these quantities. 
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Next, let us assume that at t — t, the system is allowed to abandon either 
of regimes (1. or 2.) and enter into a dissipative phase. We shall now find 
the conditions for these transitions to occur. 

Let us start with the case 1. Then, evaluating the transport equation 
immediately after the system leaves the perfect fluid regime (t ~ t + O(r)), 
we obtain that eq.(|33|) yields 



Q 



kT 



-U) 



(37) 



where (|34D and the fact that X" ~ 0(u>) ~ 0(Q), have been used. Observe 
that immediately after the system leaves the perfect fluid regime 



P r -P± = Qlj = 0{uo 2 ) 



(38) 



therefore the contribution of the shear viscosity at t — t is one order smaller 
than the pure heat conduction effects, and accordingly will not be considered 
here. 

Next, it will be more convenient to use the "conservation laws" (Bianchi 
identities) 



T.r = o 



(39) 



instead of field equations. 

The r-component of (|39|) reads 



rpr 



(p + Pr) U 

t n {l-UJ 2 ) ~ t n (l -UJ 2 ) 



t n (i - u 2 y 



t n (i 



ijj- 



+ 
+ 
+ 



4 ( ■ i n 

— \n In t H — 

t n \ t 



{p + P r )uJ Q{l+U0 2 



1 

2 



l-uo 2 l-uo 2 
(p + P r ) 2uu' 2Q'u (1 +uj 2 )2Quj'' 



rt 



2n 



(1-u 2 ) 2 

(P + P±)^ 2 



1 



uj 2 y 



Quo 



1 



which evaluated immediately after t, yields 

Q = (p + P r ) Co 



(40) 



(41) 
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Combining ( |3T| ) and (IT]), we finally obtain 



^(P + Pr) 1 



0. 



(42) 



Observe that in (pLTD and (^) we may replace P r by p, since they differ by 
terms of order 0(u 2 ). 

Now, since initially (t < t) the system was evolving as a perfect fluid 
(to — Q — 0), then a necessary condition for leaving that regime at t — t is 



kT 



r(p + p) 



(43) 



for otherwise Co = 0, and, as can be easily checked by taking consecutive time 
derivatives of ( |33| ) and (f40"|), all higher time derivatives of u; (and Q) will also 
vanish. As mentioned before, eq.(^) defines the critical point. 

Let us now consider the case 2. In this case the system (initially) is 
evolving as a quasi perfect fluid, leaving that regime at t — t. Evaluating 
( |3"3"| ) at t < t, we obtain (observe that with the choice of (^), Q > implies 
that q r points inward) 



Q = K 



V ( n 

— + Tlu{ h\nt+- ) +uoT 



t 



(44) 



then evaluating the same equation ([33]) immediately after t, and using ([44]) 

we get 

n kT - 



(45) 



as in the previous case (eq.(|37])). 

Next, evaluating (EH) just after i one obtains 



to 



(p + P r ) Q | u(p + P r ) , 4 
t n t n 



f 



or, using (|28|WBIb and <M 



4 / n \ 
+ -[h\nt+j)[(p + P r )LU-Q]=0 (46) 



Co (p + P r 



1 - 



r(p + P r ) 



-io [p + P r 



(47) 
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In the case n =constant, ( fiTD becomes (using (|28| ) and fl30|) and replacing P r 
by p) 

u = 2 -^- (48) 



1 



kT 
t(p+p) 



t 



Therefore to ^ 0, and is larger, the closer is the system to the critical point. 



4 Conclusions 

We have assumed in this paper that the zero-curvature FRW metric repre- 
sents consecutively, first a perfect (or quasi-perfect) fluid solution, and then 
a dissipative solution of Einstein equations. 

Necessary conditions for the transition from one state to another have 
been found. 

If the system is initially evolving as a pure perfect fluid, then transition 
to the dissipative regime demands that the fluid satisfies the critical point 
condition. One important question related to this case is what could be 
the physical reasons for the system to abandon the equilibrium state?. One 
possible scenario might be the appearance of a dissipative process like particle 
creation, which as is known, is formally equivalent to the introduction of 



viscous terms ||12|| . Another possibility might be the decreasing of opacity 
of the fluid from very high values, preventing the propagation of null mass 
particles as photons and neutrinos (trapping), to smaller values allowing 
for radiative heat conduction and viscosity.A similar situation occurs during 



the gravitational collapse of stars (see |13| and references therein) and the 
Kelvin-Helmholtz phase of neutron stars formation (see [T3| and references 
therein). 

Of course the full implementation of any of these, or any other, scenarios 
would require a much more elaborated setup than the one presented here. 
The use of FRW with a single dissipative fluid is obviously very limiting, 
however as we have seen, it leads to mathematically tractable equations. 

In the same order of ideas, observe that if the system is, in the initial 
phase, inflating exponentially with p + P r = 0, then as it follows from (|4l|) 
and (|45|), it will not leave that regime as long as p + P r = 0. This is an 
obvious consequence from (^)-(pTT]) which implies Q = if p + P r = 0, 



parenthetically it is worth noticing that the exit from inflation is usually 
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explained in terms of dissipative process like particle production. From the 
above, it follows that though the obtained results are independent on the 
equation of state, some specific cases are not allowed since they exclude the 
possibility of dissipation (e.g. p + P r = 0). 

If the system is allowed to evolve initially as a quasi-perfect fluid, then 
transition to dissipative regime may always occur, without requiring the crit- 
ical point to be attained, however such transition will be "faster" , the closer 
is the system to the critical point. 

Another important question to ask is if real physical systems may reach 
the critical point. Indeed, causality and stability requirements obtained from 
a linear perturbative scheme are violated, close to, but below the critical point 
(in the absence of viscosity) ||. On the other hand however, examples of flu- 
ids attaining the critical point and exhibiting reasonable physical properties, 
have been presented elsewhere [|7|,|T5|]. 

The possible solution to this apparent contradiction being the non- validity 
of the linear scheme used to obtain causality and stability conditions, close 
to the critical point (see a discussion on this point in 0). At any rate, in 
presence of both, heat conduction and viscosity, the corresponding critical 
point is beyond causality and stability conditions [0. 
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